Self-organization of atoms along a nanophotonic waveguide 
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Atoms coupled to nanophotonic interfaces represent an exciting frontier for the investigation 
of quantum light-matter interactions. While most work has considered the interaction between 
statically positioned atoms and light, here we demonstrate that a wealth of phenomena can arise from 
the self-consistent interaction between atomic internal states, optical scattering, and atomic forces. 
We consider in detail the case of atoms coupled to a one-dimensional nanophotonic waveguide, and 
show that this interplay gives rise to self-organization of atomic positions along the waveguide, which 
can be probed experimentally through distinct characteristics of the reflection and transmission 
spectra. 



In recent years, ultracold atomic gases have become rich systems for the investigation of novel many-body phenomena 
involving spatial and/or spin degrees of freedom pQ. The success of this field takes advantage of the ability to 
trap atoms in optical lattices and manipulate their interactions, enabling one to controllably emulate more complex 
condensed matter systems. In the majority of these experiments, however, the underlying lattice is static and provided 
by external lasers, which prohibits one from accessing phenomena associated with the back-action of lattice dynamics 
on spatial and internal dynamics. 

This potentially rich behavior has been explored recently using cold atoms coupled to Fabry- Perot optical cavities [2]- 
7 j. The key feature of these systems is the interplay between forces arising from photons scattered into the cavity 
by the atoms and the position-dependent atom-cavity coupling strengths. In the case of single-mode cavities, for 
example, the atomic spatial configuration undergoes a spontaneous symmetry breaking into one of two states dictated 
by the cavity boundary conditions [2j[3]. The formation of more exotic phases of matter, such as supersolids and 
quantum spin glasses, has been predicted [4j [5] and observed [6j [7] by exploiting spatial interference between cavity 
and external pump fields or multiple cavity modes. 

In this Letter, we show that atoms coupled to nanophotonic systems constitute a versatile platform for investi- 
gating similar behavior in the absence of externally imposed trapping potentials and boundary conditions set by 
cavity mirrors. This exciting frontier is motivated by recent experimental success in coupling cold atomic gases to 
nanophotonic waveguides [8-10 , and by predictions that certain atomic spatial configurations can yield remarkable 
optical response [TTTH4] , such as highly-reflecting atomic mirrors [I3j [14] . Whereas this previous work assumed the 
atoms were spatially fixed, here we show that the atoms can in fact self-organize into exotic spatial configurations 
due to the interplay of atomic internal states, positions, and photon- mediated forces [T5] . 

Our system model generally describes an ensemble of N atoms with ground and excited states \g), |e), transition 
frequency ujq, and positions zj (1 < j < A/"), which experience dipole coupling with equal strength to a single- 
mode waveguide with left- and right-propagating modes (see Fig. |TJl) . The atoms interact with each other via 
photons scattered into the waveguide modes, and effectively eliminating the associated fields produces a quantum 
spin model that describes atomic dipole-dipole interactions [14]. In particular, the reduced atomic density matrix 
obeys p = -i[H dd , p] + £ dd [p], where 

H dd = ^Sin k \Zj - Z r \cr J eg Cr J g e , (l) 

3d' 

c dd [ P ] = y, 2 °ip°\-°\°ip-p&A- (2) 
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Here, the jump operators are defined by (D± = y/Tm/4: . a J ge e Tlk ° Zj , where Tid denotes the spontaneous emission 
rate of a single, independent atom into the guided modes, fco is the wavevector at the atomic resonance frequency, 
and a 3 ge — \g 3 )( s e?\ is the lowering operator for atom j. Physically, H dd describes coherent dipole-dipole coupling 
between atoms j, j', while £dd[p] captures cooperative emission (such as super- and subradiance [16]) mediated by 
the two guided directions (I = ±). An important feature of our ID system is that the dipole interactions involve 
all pairs of atoms, and are approximately infinite in range and oscillatory in strength. This reflects the capability 
of the waveguide to channel a photon emitted by one atom without diffraction to a distant second atom, with the 
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interaction depending only on the relative phase between the atomic coherences. Corrections from causality and 
free-space propagation phase (the derivation assumes that all relevant wavevectors are close to fco) are negligible 
for realistic system sizes and atom number N [13] . However, we must account for a non- negligible spontaneous 
emission rate T' into free space, which we include via additional Liouvillian terms with independent jump operators 

6' 3 = v'T^H (i < j < N). 

We emphasize that photon-mediated forces can be extremely large in nanophotonic systems, which allows their 
effects to be prominent. Specifically, the strength Tid of the pairwise interaction in Eq. can be a significant 
fraction of the vacuum emission rate To, which corresponds to the Doppler temperature in atomic motion. As shown 
below, the long-range nature of interactions can produce forces further enhanced by a factor ~ N over a single pair. 
Current systems with nanofibers attain coupling strengths approaching Tid ~ O.lXo [8HT0], with further increases 
expected from improved waveguide design and trapping techniques. 

Now we consider the case where the atoms are driven identically by a pump field with Rabi frequency Q and 
detuning S = cj pump — c^o, as illustrated in Fig. |TJl. We focus on the semi-classical limit, in which the motion of the 
atoms is treated classically and atomic saturation is ignored. The latter assumption allows for different atoms to be 
de-correlated, e.g., (<ji<ji ) ~ (<j j )(<j j ). The evolution equations of the operator mean values satisfy 

Zj pj/m, (3) 

al e = (iS - r/2)^ e + in - ^ £ °Ce ikolZj - Z *' ] > (4) 
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(5) 



For notational simplicity, we will avoid the explicit use of brackets () to denote the mean values. Here T = + T f 
is the total spontaneous emission rate of an independent atom, m is the atomic mass, and pj is the momentum of 
atom j. An equivalent formulation of the optical forces based upon transfer matrices was derived in Ref. [15], which 
investigated the formation of ID optical lattices given asymmetric trapping fields. 

Eqs. (|3|-([5| determine the full atomic dynamics and the possibility of self-organized configurations. We now briefly 
summarize our approach to finding these solutions. A configuration is self-organizing if it is a steady-state solution to 
Eqs. ([3|-([5| and if this solution is stable with respect to small perturbations. In the "weak-scattering" limit of large 
pump detuning or low optical depth, we show that there is a unique solution that satisfies an energy minimization 
condition. Generally, however, it is not possible to formulate the problem in terms of energy minimization. Thus, we 
employ an adiabatic procedure, in which the weak-scattering solution at large detuning is continuously transformed by 
bringing the detuning toward resonance in small steps and finding the new stable configuration each time. Numerically, 
we accomplish this by introducing a small, external momentum damping rate for each atom and integrating Eqs. (J3|- 
([5| for each new detuning until the system reaches steady state. Such an external damping term is in fact needed to 
achieve true steady-state behavior, to compensate for internal anti-damping forces arising from the delayed response 
of the atomic coherences to the motion. In practice, however, we find that anti-damping can be made negligible over 
experimentally relevant time scales. 

To describe the self-organized configurations, we assume without loss of generality that Zj < Zj> for j < j' and define 
z\ = 0. While the forces in Eq. (J5| only influence the relative coordinates and the center-of-mass motion of the atoms 
in principle remains free, in practice this motion will have nearly zero velocity due to the initial cooling required to 
load atoms along the waveguide [9] [10]. It is convenient to write the positions in the form zj/Xq = rij + /j, where rij 
is an integer and fj is a fractional distance (0 < fj < 1), as shown in Fig. [IJi. Due to the periodic nature of Eqs. Q 
and ([5|, the integers rij have no observable consequence in the evolution, and the system is fully characterized by fj. 

The "weak-scattering" regime is defined as that where the fields re-scattered by the atoms into the waveg- 
uide (proportional to the last term on the right-hand side in Eq. Q) are negligible compared to the external 
driving field. A sufficient condition is A/Tid <C ^ 5 2 + (r/2) 2 , in which case the atomic coherences are equal and 
given by a J ge w ag°J = ift/(T/2 — iS). The atomic motion is then governed by a purely mechanical potential, 
#dd ~ (riD-so/2) J2j j' sin&o|£j — Zj'\, where so = |cr^| 2 - This potential is minimized when the atoms form a lattice 
with lattice constant d ws = Ao(l — ^) (see Fig. [Tj}), with corresponding energy #dd,min = —N 2 Fl ^ s ° . The N 2 
scaling reflects the strong collective forces arising from the long-range interactions within the system. The fractional 
distances take the form fj = l — (j — l)/2iV, as illustrated in Fig. [1J for the case of N = 10 atoms. 

An atomic lattice spaced by the resonant wavelength Ao = 27r//co forms a perfect mirror on resonance as N — >> oo, 
even though a single atom is mostly absorptive [14]. This phenomenon arises due to the coupling of light to a 
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superradiant spin wave S = J2j cr J ge of the atoms. As described below, for large TV the close proximity of the weak- 
scattering configuration (with lattice constant d ws = Aq(1 — ^yO) to the superradiant one gives rise to intriguing 
crossover behavior as the pump frequency is tuned closer to resonance (S — » 0) and the atomic optical depth increases. 

Starting from the weak-scattering solution at large detuning, we obtain adiabatically transformed solutions by 
changing S in small steps and integrating Eqs. ([3|-([5| at each step until the system converges to a stationary state. 
A small external damping pj = —jePj is added to the equations to facilitate convergence. Fig. |2^i depicts some 
self-organization solutions fj for a representative set of detunings, in the case of N = 150 atoms and Tid = T/4 (also 
see Fig. |2|d depicting continuous variation of 5). Compared to the weak-scattering regime, these solutions exhibit a 
variety of interesting phenomena including lattice compression, expansion, and phase slips or fragmentation. While 
we have chosen a particular parameter set for illustration here, we emphasize that our conclusions are quite general. 

We first consider the approach toward resonance (S = 0) starting from large negative detuning. Our calculations 
reveal a gradual compression of the lattice constant, d < <i ws , as S/T — » —1/2, followed by a rapid expansion 
between -1/2 < S/T < (Fig. |^) and a crossover into the superradiant regime d « Ao (nearly equal fj). A 
simple model for moderate scattering strength can be derived by noting that the atoms primarily act dispersively 
for large detunings. The phase shift imparted in transmission by a single atom can readily be calculated [HI [17], 

0t — — arctan ^f^TT^T^ ) ' wn ^ cn leads to an effective frequency-dependent wavelength A e ff « Ao(l — t /2n). For 
negative detuning, A e ff < Ao (i.e., the medium has an effective index n e ff > 1), which accounts for the smaller lattice 
constant d < d ws . Furthermore, the dispersive shift is maximized when S/T ~ —1/2. In Fig. [2^, we plot the predicted 
lattice constant <i e ff = A e ff(l — 1/2N) in the regime S < 0, as compared to the numerically obtained lattice constants. 
Our simple model reproduces well the actual behavior, even for large optical depths when one expects the system to 
be better described by band structure than an effective index [T8] . 

Approaching resonance from large positive detuning, a similar argument would predict lattice expansion, with the 
system attaining the superradiant lattice constant of d = Ao at a critical detuning of S c ~ NT m/27r. Such an argument 
is inconsistent, however, as the atoms do not behave as independent refractive index elements in the superradiant 
regime. Instead, we find that the system crosses over into a new regime around S ~ S c , where the single lattice 
fragments into two smaller lattices with a "phase slip" between them (e.g., see Fig. [2}i for detuning S = T/2). In 
particular, the left and right halves of the atomic system have approximately constant fj (indicating a lattice constant 
in those segments of ~ Ao), while a phase slip of Af w — 1/4 occurs between those two segments (corresponding to 
a 3Ao/4 separation between the two halves). Physically, each segment behaves like a collective "super-atom" with 
large dipole coupling to the guided modes. These two super-atoms are bound together by optical forces in the same 
way that two single atoms would self-organize. Specifically, the spacing between the two super- atoms is expected to 
minimize a reduced two-particle potential H^d sin ko\zR — zl\, which occurs for zr — zl = 3Ao/4 and reproduces 
the observed phase slip. The association of this behavior with superradiance can also be seen in Fig. |3^i. Here, we 
have plotted the mean atomic excited state population normalized by the independent atom result, (\(ig e \ 2 } /so, as 
a function of detuning. For positive detunings, the excited state population is strongly suppressed due to the large 
collective decay rate into the guided modes. 

Near resonance, we find that the phase slip rapidly vanishes and the two segments merge to again form a single 
lattice. Thus far, we are not able to produce a simple effective model of this phenomenon, as the strong, long-range 
interactions in this regime cause simple "mean-field" descriptions involving just a few parameters to apparently break 
down. However, we believe that the nature of this many-atom, strongly interacting system will be a rich area for 
future investigations. 

For small displacements around the self-organized solutions, the atomic motion is well-described by a set of normal 
modes (i.e., phonons). A good description of these modes can be obtained by first noting that the ratio of the atomic 
recoil frequency uo r = (hkQ)/2m to the spontaneous emission rate T is typically small. This implies that the atomic 
coherences cr ge follow the motion nearly instantaneously. Writing Eq. fi| in the form & J ge = Mjk(z(t))a ge + ifl, the 
lowest-order solution is given by 

a geMst = -iM- 1 (z(t))n. (6) 

Substituting this expression back into the force equation ([5| and linearizing around the equilibrium positions z eq 
yields a set of restoring equations pj = —Kjk(zk — z e q,k)- 

In the weak-scattering limit, Kjk is a circulant matrix, which enables the phonon spectrum to be solved exactly. 
Specifically, we find N phonon modes with wave numbers 2irj/N (0 < j < N — 1) and frequencies 



i ?r sin(7rM/") Xll/2 

2(jJ r SoTm cot ; 



27V cos(2tt j/N) - cos(tt/N) J _ 



(7) 
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The trivial j = mode denotes the free (zero-frequency) center-of-mass motion. The other modes scale with large 
TV approximately as cj p h,j ~ vWso^VTid- In Fig. pp, we plot the normalized phonon frequencies oo p h,j / V^r^o^^m 
as a function of detuning. The phonon spectrum exhibits a strong softening in the superradiant regime (Fig. |3^i). 
It should be noted that our approach to find self-organized solutions based upon adiabatic transformation does not 
preclude the possibility of other solutions. In fact, preliminary calculations suggest that the phonon softening enables 
the system to develop a multitude of other weakly stable solutions. This phenomenon will be studied more extensively 
in future work. 

Small delays in the atomic internal response compared to the motion can lead to damping or anti-damping forces. 
We can characterize this effect perturbatively by substituting Eq. (|6| into Q and iteratively finding higher-order 
corrections. In particular, we can write (J ge (t) « cr^ e ,inst + a ge,d, where the first-order correction satisfies 
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Like before, this expression can be substituted into Eq. ([5| to yield pj = —Kjk(xk—x eqi k)—LjkPk, where the matrix Ljk 
characterizes momentum damping or anti-damping. This term produces an imaginary component in the normal mode 
frequencies, 0J p h,j —> ^phj +^7ph,j 7 where 7 P h,j > (7 P h,j < 0) indicates anti-damping (damping). For sufficiently large 
atom number N, we generally find that some of the modes exhibit anti-damping. In the weak-scattering limit, the 
largest anti-damping rate behaves approximately as 7 max ~ ~ Fl ^ SoUJr . While anti-damping in principle implies that 
the system is only meta-stable absent an external cooling mechanism, in practice, the anti-damping rates can be made 
negligible compared to typical inverse trapping lifetimes for cold atoms. In the weak-scattering limit (iVTiD <C \S\), 
for example, 7 max <^ so^ r is a small fraction of the recoil frequency. In Fig. (3}d, we plot the maximum anti-damping 
rate 7 ma x versus detuning. 

Self-organization can be observed through its dramatic influence on light propagation through the waveguide. Here, 
we consider the linear reflection and transmission of an incident guided probe beam, as illustrated in Fig. [TJl. For 
concreteness, we do not account for the atomic forces or motion imparted by the probe. For example, once the atoms 
evolve to the steady-state, self-organized positions as in Figs. [2^i,b, the resulting configuration would be probed in 
a destructive fashion by a probe pulse £? pro b e (£) of duration r <C cj p h- As discussed in Refs. [HJ[T7], the reflection 
and transmission coefficients of a single atom are given by r(S p ) = — fir^ and t(S p ) = 1 + r(5 p ), where S p is the 
probe detuning. Given the coefficients of a single atom, one can efficiently calculate the reflection and transmission 
amplitudes of an array (including free-space propagation between atoms) via the transfer matrix technique [T5| fTg] . 

It is helpful to first consider light propagation through an infinite, perfect lattice with lattice constant d, in which 
case the optical modes are well-described by Bloch wavevectors [19] . Following the techniques of Ref. [18] , we find that 
the optical Bloch wavevectors q obey the dispersion relation cos qd ~ coskod — (sinkod, where ( — ~r^" l+^^/r 7 ) 2 • 
Furthermore, in analogy with purely dispersive media, we are motivated to temporarily ignore the atomic absorption 
(imaginary part of Q and look for "band gap" regions (where q is purely imaginary), in which propagation is forbidden 
due to strong interference in the multiple reflections. Defining a parameter e = 2ir(l — d/Xo) that characterizes the 
offset from the superradiant lattice constant d = Ao, we find that an optical band gap exists for probe detunings 
— < (5 gap < — Ijf^ (assuming that Vm/V ^> e). In dispersive media, a band gap gives rise to perfect reflection 
from a long lattice due to the absence of propagating modes within the system. 

Self-organization in our system into a lattice with d < Ao (e > 0) thus should manifest itself in an asymmetric 
peak in the reflectance spectrum, as the band gap occurs for negative probe detunings S p < 0. The width of this 
reflection peak should scale approximately like e _1 . In Fig.Ji^i, we have plotted the reflection spectra for the different 
self-organization configurations (as determined by the pump detuning S). Here, we have evaluated the spectra based 
on the numerical solutions of Fig. [2}d and including fully atomic absorption. We have overlaid the edges of the band 
gap — < S gap < — (dashed blue curves) in the regimes where the system is well-described by a single lattice 
constant (i.e., outside of the phase slip configuration). A clear correlation between enhanced reflectance and the 
predicted band gap is observed, with the absence of perfect reflectance attributable to the finite atom number, atomic 
absorption, and small variations of e along the atomic chain. 

In the phase slip configuration that occurs for positive pump detuning, the reflection spectra change dramatically. In 
analogy with Ref. [14] , this system can be modeled as two high-reflectivity atomic mirrors forming a high-finesse cavity 
mode, which enables its optical properties to be easily determined. Here, the spectrum becomes symmetric with respect 
to S p , with the peak reflectance and full-width given respectively by R « 1 — 4r / /(7VTiD) and #fwhm ~ ATid/a/2. 
The transition from asymmetric to symmetric spectra as the atoms cross over to the phase slip configuration is 
evident in Fig. [4^. In Fig. |4Jd, we have plotted the peak reflectance (maximized over probe detuning 5 P ) for each 
spatial configuration determined by the pump detuning S. 
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In summary, we have shown that cold atoms coupled to nanophotonic waveguides can exhibit rich self-organization 
behavior due to the interplay between atomic motion, internal states, and optical response, when the dynamics 
are treated in the semi-classical regime. However, we believe that a variety of novel quantum phenomena should 
arise as well. For example, it would be interesting to investigate how the optical fields become entangled with the 
atomic positions [20], or account for correlations that would reveal the emergence of quantum phases of matter [5]. 
Furthermore, it should be possible to tailor more exotic spin- force Hamiltonians by exploiting multi- level atomic 
structure (see, e.g., Ref. [14] for a A configuration) and by using photonic crystals [19 to engineer the dimensionality 
and geometry of the atom-nanophotonics interfaces. Beyond the exploration of many-body phenomena, we envision 
that these systems can also find applications in quantum information processing or few-photon nonlinear optical 
devices [21]. For example, it may be possible to use phonons in self-organized configurations to manipulate and 
transfer quantum information over long distances, much like in ion trap systems [22], or to exploit self-organization 
arising from an initial pump pulse to exert a nonlinear response on a subsequent probe pulse. 
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FIG. 1: a) Schematic illustration of N atoms coupled to a single-mode waveguide, taken here to be an optical nanoflber, with 
physical positions along the waveguide z\ < Z2 < ... < zn- The atoms are driven identically by a uniform pump field with Rabi 
frequency Q, and detuning 5. The atomic positions Zj/Xo = rij+fj, in units of the resonant wavelength Ao, are parameterized by 
an integer rtj and fraction < f 3 < 1 . The system is fully characterized by fj , as the integers rtj have no physical consequence 
on the dynamics due to periodic, infinite-range interactions between atoms. A weak field of detuning 5 P incident through the 
waveguide separately probes the atomic configurations, e.g., via reflectance, b) In the "weak-scattering" limit, the minimum 
energy state corresponds to an atomic lattice with lattice constant d/Xo = 1 — 1/2N and fractional distances fj = l — (j — 1)/2N, 
as illustrated versus atom index j in c) for the case of N = 10 atoms. 




FIG. 2: a) Fractional positions fj versus atom index for selected pump detunings S/T = —15, —1, —0.2,0,0.5,6. All numerical 
simulations are for N — 150 atoms and Tid = T/4. The dotted line represents the solution in the weak-scattering limit, 
b) Positions fj versus pump detuning. The vertical dashed lines denote the detunings for which the positions are plotted in 
Fig. [2^i. For clarity, only 16 representative atomic positions are shown here. Inset: zoom of the same plot near resonance, c) 
Characteristic lattice constant d/Xo versus detuning for S < 0. The circles and crosses denote the lattice constant as determined 
by the central two atoms and the average over all atoms, respectively. The solid curve depicts d e ff = A e fr(l — 1/2N), where A e ff 
accounts for the effective refractive index provided by the atomic medium. 
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FIG. 3: a) Mean excited state population (\(J ge \ 2 ) of the atoms versus dimensionless pump detuning 5/T, for the system 
considered in Fig. ^ The population is normalized by the result so = § 2 +^r/2) 2 ^ or an independently driven atom, b) Black 
crosses: phonon frequencies cj p h of the relative motional modes in the self-organization configurations, as a function of 5/F. 
Green crosses: maximum anti-damping rate 7 max of the motional modes, scaled up by a factor of 10 for better contrast. Both 
cjph and 7max are plotted in units of \^J^sqNTi^. We have used uo r = 10 _3 r for these calculations. 




FIG. 4: a) Reflectance from self-organized atomic configuration, versus pump (5) and probe (5 P ) detunings. The blue dashed 
curves indicate the edges of the band gap region £ gap , which is well-defined when a single lattice constant exists (i.e., outside of 
the phase slip configuration). The lattice constant used here is the average over all atoms. For the phase slip configuration, the 
reflectance spectrum is nearly symmetric with respect to S p , with a FWHM given by ~ NT ir>/y/2. b) Maximum reflectance 
(optimized with respect to S p ) versus pump detuning. All simulations are for N = 150 atoms and Vm/T — 0.25. 



